Non-diagonal reflection for the non-critical XXZ model 



Anastasia Doikou 



University of Bologna, Physics Department, INFN Section 
Via Irnerio 46, 40126 Bologna, Italy 



Abstract 



The most general physical boundary 5-matrix for the open XXZ spin chain in 
the non-critical regime (cosh(r/) > 1) is derived starting from the bare Bethe ansazt 
equations. The boundary 5-matrix as expected is expressed in terms of Fg-functions. 
In the isotropic limit corresponding results for the open XXX chain are also reproduced. 
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1 Introduction 



The open XXZ model is considered as one of the prototype models in describing a plethora of 
interesting boundary phenomena, and as such has attracted much attention especially after 
the derivation of the spectrum in the generic case where non-diagonal boundary magnetic 
fields are applied P-[5]. Our objective in the present study is to derive from first principles 
the most general physical boundary S'-matrix for the XXZ chain in the non-critical (massive) 
regime. 

Diagonal boundary S'-matrices for the open XXZ model in the non-critical regime were 
extracted in [6] using vertex-operator techniques, while parallel results obtained in [7] from 
the Bethe ansatz point of view (see e.g. [8l[9|,[T0]). Similarly, diagonal reflection matrices were 
derived in [HI |T2] for the critical XXZ model, corresponding to the sine-Gordon boundary 
S'-matrix for 'fixed' boundary conditions [13]. After the derivation of the exact spectrum 
and Bethe equations for the XXZ chain with non-diagonal boundaries the generic boundary 
S'-matrix for the critical XXZ chain was computed in [HI [15] , corresponding to the boundary 
S-matrix of sine-Gordon model [13] for 'free' boundary conditions. A relevant discussion on 
the generic breather boundary S'-matrix within the XXZ framework may be also found in [13] . 
Note also that analogous results regarding diagonal and non-diagonal solitonic boundary S'- 
matrices were formulated in [T6l ITT] using the so called nonlinear integral equation (NLIE) 
method [T8] . 

To extract the generic boundary S'-matrix for the non-critical XXZ model we follow the 
logic of [lM], i.e. we focus on the open chain with a trivial left boundary and a generic 
non-diagonal right boundary associated to the full i^-matrix [T3l [T9] . As also noted in 
[Ti] the main advantage of the approach adopted -considering special boundary conditions- 
is that one eventually deals with a simple set of Bethe ansatz equations similar to the 
ones of the XXZ chain with two diagonal boundaries. Thus all relevant computations are 
drastically simplified (see also [S]), and one may follow the logic described in [20l [2T| 171 [T2] 
for purely diagonal boundary magnetic fields. Ultimately, the boundary S'-matrix eigenvalues 
are extracted directly from the Bethe equations and are expressed in terms of Fg-functions 
(g = e~^) [22] as in [H [7], where only diagonal boundaries are assumed. In the isotropic 
limit g — > 1 the corresponding rational boundary S'-matrix for the XXX open chain is also 
recovered [23] . 
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2 Bethe ansatz and boundary /S-matrix 



Before we proceed with the Bethe ansatz analysis it will be useful for our purposes here to 
give the explicit expressions of the right and left boundary /^-matrices that give rise to the 
open Hamiltonian under consideration: 

N-l 

'N 



^ = E (^r^i + <i + cosh(r/) a^al^) - ^ cosh(r;) + ^^a^ 

i=l 



sinh(?7) coshfr^O ^ Ksmhir]) / , . /,\ r ■ ■ -, ^ n^ v\ . n 

+ / \ , X o'l T-r^ coshM)^? + % sinh(?76')(T^ (2.1) 

4sinh(r^0 2sinh(r70v ^' ^ ^ v / ; ly v ; 

where in the non-critical regime we are focusing here cosh(?7) > 1, also a^'^'^ are the 2x2 
Pauli matrices, and the boundary parameters ^, k, 6', are the free parameters of the generic 
i^'-matrix [l3l[T9], which will be introduced subsequently. 

To obtain such a Hamiltonian we consider the open chain constructed using Sklyanin's 
formalism [21], with left boundary oc I and right boundary associated to the general 
solution of the reflection equation [25J given in [131 [E] i-^- 

^-,y.^(^^^\^^-^'rii)\(^^^ «:e''Mn(2r7A) \ 
y Ke""^ sin(2?7A) sin[r/(A + zOle"^"^ / ' 

The latter i^-matrix has two eigenvalues given below: 

ei(A) = 2k sin[r7(A + 2p^)] sin[77(A + ip")] 

£:2(A) = 2k sin[?7(A — ^p^)] sin[77(A — «p^)] (2.3) 
where the parameters are defined as: 

— — = 2cosh[r7(p^ ±]9~)]. (2.4) 

Note that we assume here the parametrization used in [13] in the sine-Gordon context, (see 
also [14] the references therein). Such a parametrization is also quite practical within the 
Temberley-Lieb algebra framework [26]. The parameter Q appearing in (12. 2p may be removed 
by means of a simple gauge transformation, that leaves the XXZ i?-matrix invariant, and 
henceforth we consider it for simplicity to be zero (see also [T3l [T4]). The matrix ( 12. 2p 
may be easily diagonalized by virtue of a constant (A-independent) gauge transformation: 

diag(£i(A), £2(A)) = M-\v'. p-) K{\) Mip+, p~) (2.5) 

where M is defined as: 
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Note that the above transformation modifies dramatically the XXZ i?-matrix, so it is not 
possible to simply implement a global gauge tranformation changing the basis in order to 
diagonalize the open transfer matrix as in e.g. [27]. It is also worth pointing out the 
similarity between the matrix A4{p'^, p~) and the local gauge transformation employed for 
the diagonalization of the open XXZ transfer matrix with non-diagonal boundaries [T[|28[|29]. 

We recall now the exact Bethe ansatz equations for the open XXZ chain in the case of 
a right non-diagonal boundary and a left trivial diagonal. The Bethe equations in this case 
reduce to the following simple form (see also [Ti]): 

sm[rj{X, - f (2p+ + 1))] sm[rj{X, - i(2p- + 1))] cos[77 (A, + f )] / sin[r7(A, + |)] y^+^ 
sin[r; {Xi + f (2p+ + 1))] sin[?7 (A, + f (2p- + 1))] cos[7] (A, - f)] \sm[r] (A, - f)]/ 

= _Y\ ^^^[^(^^ ~ + ^)] sin[?7(A^ + Xj + i)] 
jj^ sin[77(Ai - Xj - i)] sin[r7(Ai + Xj - i)]' 

We consider here, without loss of generahty i] > 0, > | and Re (Aq) £ [0 , ^] Aq, 7^ 0, ^ 
(see e.g. [20] for details on this restriction). For relevant results on various representations 
of Ugisk) see [281E91I30]. 

As pointed out in [H] the integer M is associated to a non-local conserved quantity S, 
which has the same spectrum as 5"^ (for more details we refer the interested reader to [TH [28] 
and references therein) i.e. 

N , , 

M=--5„ (2.8) 

the subscript e stands for the eigenvalue. 

Our objective now is to explicitly derive the physical boundary S'-matrix, and in par- 
ticular the relevant overall physical factor, which provides in general significant information 
on the existence of boundary bound states. We define the boundary 5'-matrices by the 
quantization condition [TUl 120] 

gi2p(A)JVj^+ j^- _ {\ _ ^2.9) 



A is the rapidity of the 'hole' -particle-like excitation, and p(A) is the momentum of the hole. 

The density of a state is obtained in a standard way from the Bethe ansatz equations 
after taking the log and the derivative ^ [TOl [20l [lH [7] . More precisely the Fourier transform 
of the density for the one-hole state turns out to be: 

A- (, ,\ 0?(. ,\ I Q-2(^) /Juj\ I „-iuj\\ 

as[uj) = 2e{u) + —-—-—{e +e ) 
Jy 1 + a2{uj) 



+ 



1 



Nl + a2iuj) 



(2.10) 
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where we define the following Fourier transforms 

a„(u;) = e-^''"l-l, 6„M = (-)- a„H, 



ai[uj) 



1 



1 + 02(0;) 2cosh(f) 



(2.11) 



e(A) corresponds also to the energy of the particle-like excitation. The similarity of the latter 
formula (12.101) with the one obtained in the case of two diagonal boundaries [20l [7j is indeed 
noticeable. This is a crucial point enabling a simplified derivation of the boundary .S-matrix. 
In our case both terms depending on p"^ are assigned to the right boundary, otherwise one 
follows the logic of the fully diagonal case (see e.g. [20l ITj). 

The boundary matrix K^, of the generic form fl2.2p . has two eigenvalues ki 2, whereas 
the left boundary matrix is trivial K+(A) = A;o(A)I. From the density (12.101) . the quantiza- 
tion condition (12. 9p . and recalling that e(A) = 2^^^^ we can explicitly derive the quantities 
^0, ki^2 (see e.g. [Sni EH [7j for more details). Actually, the eigenvalues ki(A, p^) and 
k2(A, p^) may be seen as the boundary scattering amplitudes for the one particle-like exci- 
tation with S = +1 and S = — |, respectively (see also [H] for more details) 

We first compute the eigenvalue ki, which is expressed in terms of the rg(a;)-function, 
-the g-analogue of the Euler gamma function- {q = e"'') defined [22] as 

:i - q'^') 



r,(x) = (i-?)i-^n 

j=0 



< g < 1. 



(2.12) 



Using also the g-analogue of the duplication formula [22] 

,1 
'2 



r,(2x) r,.(-) 



+ 



2x-l 



T„2{x) r„2(a;+ -), 



(2.13) 



we obtain the following result for the first eigenvalue ki(A, p^ , p ) (up to a constant phase 
factor) : 



ki(A, p^, p ) = 2k, sin r] + -{2p'^ — 1] 
X A;o(A) /ci(A, p+) fci(A, p'] 

where we define: 



sm 



r]{X + -{2p- 



ko{X) = q 



-4iA 



( ^ + i ) r^s ( f + 1 



Tos ( Y + I 



(2.14) 



(2.15) 



ki{X, x) 



{2k) 



+ i(2x - 1)) r,4 (f + i(2x + 1)) 



sm 



77 A - 1(22;- 1) 



r,4 f + i(2x-i)) r,4 ^ + i(2x + i) 



(2.16) 
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We turn now to the computation of the second eigenvalue k2(A, p^), corresponding to a 
one-hole state with 5 = — |. We implement the 'duahty' transformation [211 E]) which 



modifies the boundary parameters — > —p^ in the Bethe ansatz equations (12.71) . This 
transformation is the equivalent of deriving the Bethe ansatz equations starting from the 
second reference state (the analogue of the 'spin down' state) [301 |3T]. Then we conclude for 
the second eigenvalue: 



ki(A, p+, p- 
k2(A, p+, p- 



sin 


r;(A+f(2p+-l)] 




sin 


r^(x+i{2p--l)) 




sin 


"r/(A-f(2p+-l)) 




sin 


"r^(A-f(2p--l)] 





(2.17) 



An alternative way to extract the second eigenvalue is instead of the 'kink' state with 5 = — |, 
-after implementing the duality transformation- to consider the anti-kink state consisting 
of a hole and a two-string state. Such configurations have been utilized in deriving the 
kink-antikink scattering amplitudes in the bulk XXZ model (see e.g. [7|) as well as in open 
XXZ chain with the most general boundary conditions [15], where the 'duality' p^ — > — 
cannot be implemented for the derivation of the second eigenvalue of the boundary S-matrix. 
Notice that the term depending on the boundary parameters (12.161) is 'double' compared to 
the diagonal case studied in [6l [7] . Analogous phenomenon occurs in the open critical XXZ 
chain [TUfTS] and the sine-Gordon model [13]. It is straightforward to see that in the diagonal 
limit we recover the results of [HI E]. Also, in the isotropic limit q 1, Tq{x) ^{x) and 
the trigonometric functions turn to rational, hence the generic rational reflection matrix for 
the open XXX spin chain is easily recovered (see also [23]). 

It is finally convenient to rewrite the two eigenvalues in terms of 'renormalized' boundary 
parameters defined as: 



p 



± 



1 
2 



mod( 



(2.18) 



then the similarity between (12.171) and the ratio of the 'bare' eigenvalues (12. 3p becomes appar- 
ent. We have actually derived the physical boundary iS- matrix up to a gauge transformation; 
indeed the S'-matrix of the generic form (12. 2p may be reproduced by: 



K(A, p+, p-)=M{p+, p~) diag(ki(A), k2(A)) M-\p+, p 



(2.19) 



A4 is defined in (12.51) . This concludes our derivation of the general boundary S'-matrix for 
the open XXZ chain in the non-critical regime. 
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